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The tensor renormalization group study of the general spin-S Blume-Capel model 
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We focus on the special parameter point D = 2J of the general spin-S Blume-Capel model on 
the square lattice. The phase transition behaviors with the temperature variation are discussed by 
the newly developed tensor renormalization group method. For the case of the integer spin-S, the 
system will undergo 5* first-order phase transitions with the successive spontaneous breaking of the 
magnetization M = S, S ~ 1, ...0. For the half-integer spin-S, there are similar S — 1/2 first order 
phase transition with M = S,S — l,...l/2 stepwise structure, in addition, there is a continuous 
phase transition due to the spin-fiip Z2 symmetry breaking. For the first order phase transitions, 
the location of the critical temperature is the same, independent of the value of the spin-S. It is a 
different picture compared to the previous results Tc = G. 



PACS numbers: 05.10.Cc, 05.50.-|-q, 75.10.Hk 

The Blume-Capel modeli and the extended Blume- 
Emergy- Griffiths model^ have attracted the general 
interests in the several decades. The Monte- 
Carlo algorithni^ii, finite-size scaling^ and mean- field 
approximation^ demonstrate the rich phase diagrams 
of this model. The richness originates, not only from the 
spin-flip Z2 symmetry breaking, but also from the den- 
sity fluctuations with S + l(for integer 5) or 5' -I- 1/2 (for 
half-integer S) optional values^ for Sf . Furthermore, 
It is a characterization for the experimental 3He-4He 
mixtures®'^ and metamagnets^*', which inspired the vec- 
tor version of this modeliiii^. 

The Hamiltonian of the Blume-Capel modelA with gen- 
eral spin-S is 



H = -J J2 S,S,+Dj2sf, 



(1) 



<4J> 



where the spin variable S takes 2S +1 values: (—5, —5* + 
1, S* — 1, S*). J is the coupling constant and D is the 
strength of the single-ion anisotropy. The sum of the first 
term runs over all the nearest neighbors. When D = 0, 
this model is reduced to the classical Ising model. When 
D goes to the positive infinity, the energy favorable state 
is the one with the full occupation of the smallest spin 
component. For the integer spin cases, the component 
is S* = 0. For the half-integer cases, the component is 
S = ±1/2. 

If we look at the hamiltonian of any bond linking two 
sites i, J, then we have 



H,j=^JS,Sj+D{Sf + S])/q 



(2) 



Here, q is the coordination number depending on the lat- 
tice structure. This formula can be rewritten as 

i/,,, = J{S, - S,f/2 + {D/q - J/2){Sf + (3) 

from which, it leads to the following conclusions for 
the ferromagnetic coupling(J > 0): when D > qJ/2, 



Si = Sj ~ l/2(— l/2)(depending on the spontaneous 
breaking), respectively for the integer and half-integer 
spin-S cases; when D < qJ/2, the configuration of the 
ground state is Si = Si = max(S') for any spin-S cases. 
As a result, we have a special situation, i.e., D = qJ/2, 
where the ground state is Si — Sj with (25* -I- l)-fold 
degeneracy. 

Without loss of the generality, we focus on the square 
lattice hereafter. Then, q — -i and D = 2J is the 
special parameter. The phase boundaries of the square 
lattice2.-^il. end at (T/J,D/J) = (0,2), i.e., for D = 2 J, 
the critical temperature is Tc = for the general spin 
cases. 

Motivated by the special parameter point and the 
thought how the phase boundary approaches the ending 
point, we calculate the temperature dependence behav- 
iors of this model in the case of I? = 2 J by the recently 
developed tensor renormalization group algorithm For 
a classical lattice model with the local interactions, 
the partition function can be represented as the tensor 
productii^. 



Z = Tr[]r,,,^ 



Uidi 



(4) 



the configuration of the ground state is Si ~ Sj 



0, 



where i runs over all the lattice sites and Tr is to sum over 
all bond indices. The local tensor T is defined at each 
lattice site as shown in Fig. [U and the indices l,r,u,d 
mean the left, right, up and down directions respectively. 
Here, the initial tensor T is a order-4 tensor due to the 
four indices. The initial dimension do = (25' -|- 1) of each 
order is the degree of freedom for the spin. The detailed 
information about the construction of the tensor can be 
found in Ref.^^. 

To contract the tensor network, i.e., trace over all sites, 
we will face the exponential increasing problem of the 
dimension of each order of the tensor. After n-th con- 
traction along the horizontal and vertical directions, the 
dimension varies by dn = d^^i- The contraction process 
will make the dimension go beyond the calculability. In 
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FIG. 1: (Color Online)A graphical representation of the ten- 
sor network in the left panel; the four-order tensor defined at 
the site in the right panel. 



2007, Levin and Navei^ proposed a cutoff scheme for dn 
by the singular value decomposition(SVD), to get a good 
approximation for the partition function calculation, by 
which the thermodynamic properties of the system can 
be obtained consequently. 

Recently, we proposed an coarse graining tensor net- 
work renormalization group(RG) based on the high order 
singular value decomposition(HOSVD)i^, which provides 
an accurate but low computational cost technique for 
studying two- or three-dimensional (3D) lattice models. 
The following results are all from the newly developed 
KG scheme. Hereafter, the coupling constant J is used 
as the energy unit, and Kb is set as 1. 




FIG. 2: (Color Online)The temperature dependence of 
the magnetization(black square) and occupation number of 
N{S = 0)(red circle), here, S = 1,D = 2J. 



like continuous phase transitions from the ferromagnetic 
to paramagnetic phase. With the temperature increas- 
ing, the components N{S = 1, —1) come back to the 
system again with the equal-weight occupation, which 
means that the system remains in the paramagnetic state 
of M = 0. In the limit of the high temperature, the 
occupation numbers of the three different components 
S = 1,0,-1 all go to 1/3. 



I. THE CASE OF 5 = 1 

In the case of S* = 1, there are three optional values: 
±1,0 for each spin variable. It is the case catching the 
most attention. On one side, there exists the tricritical 
point enriching the phase diagram. On the other side, it 
is affordable computational cost for the different numer- 
ical methods. 

In the limit of infinite D, the energy favorable state is 
N{S = 0) = 1, which means that all sites are occupied 
by 5* = 0. We can name S' = as the vacancy or hole, 
the state full of holes is just the hole condensed phase, 
whose quantum correspondence is the one-dimensional 
5 = 1 quantum model with single-ion anisotropy of our 
previous paper—. The hole condensed states also emerge 
in the other general integer spin-S cases. 

As is shown in Fig. [21 the system undergoes the first or- 
der phase transitions accompanying the emergence of the 
hole condense phase as the temperature varies. The mag- 
netization jumps from 1 to 0, whose location is consis- 
tent with the kink of the occupation number of the holes 
N{S = 0). The critical point reads Tc = 0.123. For the 
terrace of the magnetization M = 0, the starting point 
is the condensation of the holes, which determines the 
first order phase transition behaviors other than Ising- 



II. THE CASE OF S = 2 

Here, there are five optional values: 0,±1,±2 for the 
spin variable. Compared to the case of S* = 1, there is 
one more optional pair S — ±2, which renders one more 
jump of the magnetization from AI — 2 to M — 1. The 
system undergoes the first order phase transitions twice 
with the temperature increasing. 

The first platform M = 2 corresponds to the full occu- 
pation of N{S = 2) = 1 after the spontaneous breaking 
of the pair S = ±2. The second platform M = 1 emerges 
with the disappearance of N{S — 2, —2) and full occupa- 
tion of N{S = 1) = 1 after the spontaneous breaking of 
the pair S — ±1. The second jump of the magnetization 
comes with the hole condensed phase N{S = 0) = 0. 
With the temperature further increasing, the occupa- 
tion of the hole drops down continuously. The pairs 
S = ±1 and S = ±2 come back to the system succes- 
sively. N{S = 0) + N{S = 1,-1) ~ 1 keeps until the 
considerable temperature. Accompanying the successive 
disappearance of S" = ±2 and S = ±1, the two criti- 
cal points are located at Tc^ — 0.123 and Tc^ = 0.365 
respectively. 
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FIG. 3: (Color Online)The magnetization(black square) and 
occupation number of A''(5' = 1, — l)(red circle) as a function 
of the reduced temperature T/J, here, S = 2, D = 2 J. 
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FIG. 4: (Color Online)The temperature dependence of the 
magnetization(black square) and occupation number A''(5' = 
1/2, -1/2) (red circle), here, S = 3/2, D = 2J 



III. THE CASE OF 5 = 3/2 

For S = 3/2, there are four optional values : 
±3/2, ±1/2 for each spin variable. Similar to the case 
of 5 = 1, the system exhibits the spontaneous breaking 
which start from the magnetization M — 3/2. The first 
jump of the magnetization from M = 3/2 to M = 1/2 
results from the competition of the occupation of two 
pairs between S = ±3/2 and S = ±1/2. The sponta- 
neous breaking brings us the full occupation change from 
N{S = 3/2) to N{S = 1/2). The discontinuity of the 
physical quantities mean the first order phase transition. 

When the temperature increases further, the Ising-like 
continuous phase transition shows up due to spin-flip Z2 
symmetry breaking. The results support the physical pic- 
ture of the separated first order and second order phase 
transition line in Refj^. By the behaviors of the magneti- 
zation, the two critical points are located at Td — 0.123 
and Tc2 = 0.649 respectively. The location of Tc2 is close 
to the data shown in Fig. 2 in Ref.''. 



IV. THE CASE OF S = 5/2 

There are six optional values: ±1/2, ±3/2, ±5/2 for 
each spin variable. 

When the temperature increases slowly for the zero 
point, the system spontaneously breaks into the state 
M = 5/2, with the full occupation N{S = 5/2) = 1. 
With the temperature further rising, the system en- 
ters into the state M = 3/2 with the full occupation 
N{S = 3/2) = 1, then becomes the state M = 1/2 with 
the full occupation N{S ~ 1/2) — 1. The stepwise break- 
ing is clearly shown in the profile of the occupation num- 
ber N{S = 3/2,-3/2). At the first terrace M = 5/2, 
the degrees of freedom S = (3/2, -3/2), (1/2, -1/2) are 
gone. On the way to increase the temperature, the first 
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FIG. 5: (Color Online)The temperature dependence of the 
magnetization(black square) and occupation number A''(5' = 
3/2, -3/2)(red circle), here, S = 5/2, D = 2J. 



jump of the magnetization M and the occupation number 
N demonstrates the disappearance of N{S = 5/2, —5/2) 
and full occupation of N{S = 3/2, —3/2). The sponta- 
neously breaking of the spin pair S = (3/2, —3/2) brings 
N{S = 3/2) = 1. Then the second jump of M and N rep- 
resents the ensuing replacement of = (S* = 3/2, —3/2) 
hy N — {S ^ 1/2, —1/2). The spontaneously breaking of 
the spin pair S = (1/2, -1/2) brings N{S = 1/2) = 1. 

However, with the temperature further increasing, M 
drops to continuously, accompanying the continuous 
phase transition, and the system enters into the para- 
magnetic phase. Here, the occupation number N(S — 
3/2, —3/2) increases continuously from 0. A further cal- 
culation shows that N{S = 1/2) = N{S = -1/2), N{S = 
3/2) = N{S = —3/2) in the paramagnetic phase, where 
N{S = 1/2, -1/2) ± N{S = 3/2, -3/2) =i 1 keeps until 
the considerable temperature until the equal-weight dis- 
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tribution of all the degree of freedom in the high tem- 
perature limit. The three critical points read T^^ — 



0.123, = 0.366, Tc. 



0.685. 



Reviewing the data about the critical points, we find 
that the critical temperature for the first phase transition 
of first order is fixed at Tci = 0.123 for the above four 
cases 5" = 1,3/2, 2, 5/2, and the second one of first order 



is = 0.365 for S 



2 and = 0.366 for S = 5/2 



respectively. The location of the first order phase transi- 
tion is independent of the value of the spin-S. It is wield 
at first sight, why does it happen? 

Let us turn to the principle about the minimal free 
energy, F = U — TS. The phase transition lies in the 
competition between the internal energy and the en- 
tropy. The density fiuctuations of the optional values 
bring the change of the internal energy at the interface 
of the two different clusters composed of single spin value. 
It changes the possible configurations of the occupation 
simultaneously, which is exactly the entropy. The first or- 
der phase transition occurs with the replacement of the 
optional spin values related to the successive spontaneous 
breaking, which is irrelevant for the value of spin value S. 
The interval of S is always 1 no matter that it is the inte- 
ger or half-integer case, which induces the same internal 
energy difference. The balance of the free energy is the 
same for the different spin-S cases. For these low tem- 
perature transition points, the relative Boltzmann weight 
g-^AB .j-j^g variation of the spin value is quite 

small. As a consequence, we cann't see the appreciable 
thermal fiuctuations. 

As was pointed out in Ref.^-, a symmetry breaking 
phase with n degenerate states is represented by fix-point 
tensors, which is a direct sum of n dimension-one triv- 
ial tensor. XI is a visualized parameter of the fix-point 
tensor with the following definition. 



XI = 



Idrd 



(5) 



and its graphical demonstration was referred to Fig. (13) 
in Ref JI. 

For the integer case of S" = 2, the visualization XI 
bears the step-wise structure: 5,3,1. The 5-fold degen- 
eracy associated to T^^ T^^ rpTRi jg represented by 
the first platform valued 5. T^^ corresponds to the spin 
fiip of the pair (s, —s) and T^^^ comes frome the degree 
of freedom 5 = 0. With the temperature increasing, the 
system undergoes two phase transitions. The first jump 
of XI from 5 to 3 and T^^ originating from S — ±2 is 
gone, where N{S = 1) = 0. The platform of XI = 3 cor- 
responds to T^^^T"^^^ referring to = 0, 1, -1. The 
platform of XI = 1 is the consequent vanishing of 
from another pair S = ±1, where N{S = 0) = 1. 

For the half-integer case of S — 5/2, the visualization 
XI bears the step- wise structure: 6, 4, 2, 1. The 6-fold de- 
generacy comes from the three pairs ±5/2, ±3/2, ±1/2, 
corresponding to T^^ T^^ T^^ . The first two jumps 
of XI represent the two first order phase transitions, 
which are similar to the case of 5* = 2 accompanying with 




FIG. 6: (Color Online) The step- wise structure of the visual- 
ization XI with three platforms (5, 3, 1) as a function of the 
reduced temperature T/J. here, S — 2, D = 2 J. 
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FIG. 7; (Color Online)The step-wise structure of the visual- 
ization XI with four platforms (6, 4, 2, 1) as a function of the 
reduced temperature T/J. here, S — 5/2, D = 2J. 



the successive vanishing of the occupation pairs. The dif- 
ference lies in the last jump of XI from 2 to 1, which is the 
continuous phase transition with Z2 spin-flip symmetry 
breaking. The gap of XI is 2 for all the first order phase 
transitions and 1 for the last continuous phase transition 
respectively. 

We can come to the general conclusion that the vi- 
sualization parameter XI bears the step-wise structure: 
25"+ 1, 25* -1,..., 3,1 for the integer spin S, 2S+1,2S- 
1,...,2,1 for the half-integer spin S. It is tricky to lo- 
cate the position of the critical point by the step of XI 
because of the subtlety of XI due to the cutoff in the 
contraction process, however, the integer platform of XI 
associated to the degeneracy of the system is intrinsic, 
which comes from the deep insight that the fix-point of 
the tensor representation corresponds to the fix point of 
RG flow associated to the phase transition. It does open 
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the new window to observe the phase transition. The im- 
provement and exploration about the visuahzation like 
XI from the tensor itself is still on the way. 

In summary, we discuss the phase transitions of the 
Blume-capel model on the square lattice using the re- 
cently developed tensor network renormalization group 
scheme. We focus on the special situation with the high 
degeneracy, D = 2 J, where the bond Hamiltonian is the 
form of the square sum. The calculation shows, for the 
first order phase transition, the different phase boundary 
ending points successively locate at Tc„(n=l,2,...) on 
the vertical line D = 2 J. Here, max(n) = S for integer 
spin-S and max(n) = 5 — 1/2 for half-integer spin-S, Un- 
like the results that the ending points always touch at 
{T/J,D/J) = (0,2) from the mean-field approximation 
and finite-size scaling method. We fix the first two crit- 



ical points = 0.123, = 0.36(5) of the first order 
phase transitions for the general spin-S cases. Through 
the magnetization and the occupation number of differ- 
ent optional spin values, the first order phase transition 
behaviors are clarified, accompanying the successive dis- 
appearance of the optional spin pairs. The visualization 
parameter XI associated to the degenerate states pro- 
vides the further evidence for the phase transition pro- 
cess. Compared to the integer spin-S cases, there is one 
more Ising-like continuous phase transition with spin-flip 
Z2 symmetry breaking for half-integer spin S-cases. 
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